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Abstract 

We investigate SU (2) Skyrmions in hyperbolic space, by computing numerical 
solutions of the nonlinear field equation. We first demonstrate the link between 
increasing curvature and the accuracy of the rational map approximation to the 
minimal energy static solutions. We investigate the link between Skyrmions with 
massive pions in Euclidean space and the massless case in hyperbolic space, by 
relating curvature to the pion mass. Crystal chunks are found to be the minimal 
energy solution for increased curvature as well as increased mass of the model. 
The dynamics of the hyperbolic model are also simulated, with the similarities and 
differences to the Euclidean model noted. 
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1 Introduction 


The Skyrme model [T] is a (3+l)-dimensional nonlinear theory of pions that admits 
topological soliton solutions, called Skyrmions, which represent baryons. This has been 
well studied [2] with solutions calculated for a large range of topological charges [3] . The 
addition of a mass term has little effect on solntions of low baryon nnmber, which continue 
to form shell like structures. However for larger charge solutions, a mass term starts to 
favour minimal energy solutions formed of finite chnnks of a Skyrme crystal ilElle]. 

It has been demonstrated that there is a surprising similarity between Skyrmions with 
massive pions in Euclidean space and the massless case in hyperbolic space [7]. The 
cited paper also outlines a method for constructing Skyrmions with massive pions from 
instanton holonomies, by first modelling a hyperbolic Skyrmion by taking holonomies 
along particular circles in [H] and applying a mapping relating hyperbolic cnrvature 
and Euclidean mass to produce the Euclidean Skyrmion [Zj. 

This posits that there could be a geometrical underpinning to the standard mass term, 
traditionally nsed in the Skyrme model. This suggests that understanding Skyrmions in 
hyperbolic space and the affect that cnrvature has, may shed some light on Skyrmions 
with massive pions in Enclidean space. Most notably there are certain properties for 
Skyrmion solutions in Euclidean space, that only occur once the mass term is turned on, 
or exceeds a certain threshold. Namely, the formation of crystal chunk solutions, as the 
global minima, for higher charge systems that exceed the threshold mass. If some similar 
behaviour were to be observed for massless solutions in hyperbolic space, it wonld snpport 
this geometric link. In fact, it will be demonstrated that the map linking the cnrvature 
of hyperbolic Skyrmions with massive Euclidean solutions, can be used to predict the 
global minima solution. 

We will also examine the dynamics of Skyrmions in hyperbolic space, demonstrating 
that they scatter along geodesics, with maximally attractive channels corresponding to a 
relative rotation through an angle tt, about an axis orthogonal to the connecting geodesic. 


2 The Model 

The Lagrangian density for an SU(2) valued Skyrme field U{t,x) is given by, 

jC = ArriR.Rn + ^Tri[R„ R,][R^, - mlTr{U - I 2 ) (2.1) 

where Ri = {diU)W is the right su{2) valned cnrrent. The associated energy for a 
static Skyrme field U{x) defined on a general Riemannian manifold M with metric = 
Qijdx^dx^ is 

E = [-\Tr{RiR^) - ^Tr{[R,,R,][R\R^]) +m^Tr{l - U)^^d?x (2.2) 
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where g is the determinant of the metric. Note that both of the above expressions have 
the parameters preceding the hrst two terms scaled out. m is the tree-level mass of the 
pions, which can be observed by using the SU{2) nature of the field and writing the 
equation in terms of pion fields U = a + iiv ■ t, where r is the triplet of Pauli matrices 
and TT = (vTi, 7r2, TTs) the triplet of pion fields. 

Much work has been done on the solutions to this equation for Euclidean space M = 
upto topological charge 108 mu- However we are interested in considering Skyrmion 
solutions in hyperbolic 3-space M = which is the space with constant negative 
curvature —The metric of takes the form, 

ds^ (H^) = dp^ + ^ ^2.3) 

where p is the hyperbolic radius. If we take the limit of zero curvature, we recover the 
Euclidean metric, with the hyperbolic radius equal to the standard Euclidean radius p = 
r. We will also make use of the standard Poincare ball model for displaying results. This 
can be obtained from the above metric by a simple radial transformation p = ^ 

to give the following metric. 


ds\K) 


(1 - 


(2.4) 


Hence our space can be modelled by a sphere with a boundary at infinite hyperbolic 
radius given by i? = ^ (though our plots will always be scaled to an equivalent size). 

The vacuum for the massless theory is any constant U, however the inclusion of the mass 
term m > 0 gives the unique vacuum to be 17 = I 2 . We will impose the boundary 
condition U —>-12 as p—>- 00 , which is required for finite energy. This gives us a map 
U : U { 00 } = S^, and hence a topological charge as an element of the 3rd 

homotopy group, equivalent to an integer B G vr3(5'^) = Z, 


B = - 


2471 ^ 


tijkTr{RiRjRk)Sx. 


(2.5) 


3 Approximations 

There are a few approximations for Skyrmions with massless and massive pions. The 
rational map approach will be the most useful in this paper. The angular dependence 
of the solution is approximated to be a rational map between Riemann spheres [9]. On 
extension to massive pion solutions, it is found that only shell-like approximations can be 
closely approximated. While multi-shell like solutions have been modelled in an attempt 
to form more crystal like solutions (III. they are poor approximations to the full mini¬ 
mal energy solutions. They can be useful for initial conditions in numerical simulations 
however. 
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(a) energy density plot 




(b) profile function f{p) (c) energy for increasing curvature 
(k = l,m = 0) (m = 0) 


Figure 1: B = 1 static hedgehog solution, (a) energy density plot in Poincare ball, where 
grey shaded region represents the boundary of hyperbolic space, (b) profile function f(p) 
for K = 1, m = 0, (c) energy for increasing curvature, for m = 0. 


3.1 B=1 

In Euclidean and hyperbolic space the single Skyrmion solution can be reduced to solving 
an ODE, using the hedgehog ansatz. This is known as a hedgehog solution due to its 
radial nature, as can be seen in figure The field is given to be 


U = exp{if{p)x ■ r), 


(3.1) 


where x = (sin d cos 0, sin d sin 0, cos 0) is the unit vector in Cartesian coordinates, f(p) 
is a monotonically decreasing radial profile function with boundary conditions /(O) = tt 
and /(oo) = 0. Substituting this into the energy in (2.2) we get a radial energy of the 
form. 


E = 


1 

Stt 


/ 


/2 


sinh^ Kp 




+ 2(r + l)sinV + - 


2 • 4 r 

K Sin j 


sinh Kp 


+ 2m' 


, sinh^ Kp 


K^ 


(1 - cos/) 


dp 

(3,2) 


The profile function /(p) can then be found by minimising the above energy and is also 
shown in figure[l]for k = l,m = 0. This yields a function with an exponential asymptotic 
decay for m = 0, 


/ ~ (3.3) 

This takes a similar form to that of massive Euclidean Skyrmions (k = 0) / ~ but 

dependent on the curvature rather than the mass of the theory. This suggest a relation 
between curvature and mass. In fact it is found that if you select the correct curvature, 
you can produce an extremely similar profile function for any Skyrmion with massive 
pions in Euclidean space. See [7j to observe the graph showing the relation between k 
and m. 
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3.2 Shell-like multisolitons 


Shell-like solutions can be well approximated by the rational map ansatz. In hyperbolic 
space this takes the following form, 


U (p, z) = exp 


if{p) 

1 + \R\^ 


2R \ 
2R ) 


(3.4) 


where z = e*‘^tan(|) is the Riemann sphere coordinate and R{z) is a degree B rational 
map between Riemann spheres. Substituting this ansatz into ( | 2 . 2 [ ) we get the following 
radial energy, 


E= — 
Stt 


/ 


/2 


sinh^(Kp) 




+ 2B{f^ + l)sm^ f + 1- 


sin^ / 


sinm(Kp) 


-|- 2m' 


,sinh^(Kp) 




(1 - cos f)jdp, 
(3.5) 


where 


X = 


dyr 




12 

dR 

\ 4 2idzdz 

Vi + |fip 

dz 

J (1 + 


2)2 


(3.6) 


X is an integral to be minimised by the choice of rational map R{z). Note that X is 
independent of k and hence the values match those in Euclidean space. The minimal 
values of X and the associated rational maps can be found in [3] for a range of values of 
B. Note that the earlier hedgehog ansatz is recovered for B = 1, where R = z is the 
minimising map, with X = 1 and (3.5) rednces to (3.2). 


This approximation will be used in various way to form initial conditions for the numerical 
computations presented later. We will also investigate how curvature affects the accuracy 
of the approximation. 


4 Static Solutions 


4.1 Shell-like Static Solutions 


The static eqnations that follow from the variation of ( 2 . 2 ) were solved using a time depen¬ 


dent 4th-order Runga-Kutta method to evolve the time-dependant equations of motion 


that follow from the relativistic lagrangian ( 2 . 1 ), cntting the kinetic energy whenever 


the potential increased. The grid was modelled using the Poincare ball model of radius 
on a cubic grid with (201)^ grid points and lattice spacing (for the standard k = 1 ) 


K 


-1 


Ax = 0.005. Spatial derivatives have been approximated using a 4th-order finite differ¬ 
ence method. We mnst fix the bonndary at i? = k,~^ to be the vacnnm at spatial infinity 
I 2 , to ensnre finite energy. For all onr simnlations the topological charge, when 


Noo = 


5 



















Table 1: The energy for soliton solutions (E) and rational map ansatz (Er) with k = 
l,m = 0 


B E E/B Er Er/B % difference figure 



computed numerically, gives an integer value to five significant figures, indicating the 
accuracy of the results. 


Two forms of initial condition were considered. The rational map ansatz shown in (3.4) 
and the product ansatz U{x) = Ui{x)U 2 {x), which was used to place lower charge solitons 
at various well separated positions about the grid. 


The first eight shell-like static solutions for k = l,m = 0 can be seen in figureThese 
solutions take a similar form to the Euclidean solutions of the same charge, with a few 
subtle differences. The faces of the polyhedron now appear to take the form of geodesic 
surfaces (a surface that contains curves belonging to the set of geodesics within the global 
space). Additionally, translating the solutions about the grid alters the apparent shape 
and means that fines of symmetry fall along geodesics of the space. This can be seen in 
more detail in the analysis of the B = 8 solution in figure The crystal chunk solution 
clearly demonstrates a bowing of the line connecting the two B = 4 solitons, this line is 
found to be a geodesic of the space. 


If we look at the energies displayed in table [T] we can see the expected trend in energies 
for increasing charge. We also observe how the energy of a given charge solution scales 
with curvature in figure 

We now compare the approximation from the rational map ansatz to the minimal energy 
solution for topological charges 5 = 1 to 8. The results for 5 = 2 can be observed in 
figure 1^ We note that the rational map gives a very good approximation up to 5 = 4. 
The fraction Er/E, where Er is the energy of the rational map approximation and E 
is the full numerical minimal energy, seems to stay relatively constant throughout an 
increase in curvature. We can’t say if this trend will definitely continue, however if it 
does, then the rational maps will remain a good approximation for all values of curvature, 
as long as the solutions are shell-like, but the rational map approximation breaks down 
if the solutions begin to become non shell-like. 
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(a) S = 1 


(b) S = 2 


{c) B = 3 





{d) B = 4 


(e) B = 5 


(f) B = 6 




ig)B = 7 (h)B = 8 


Figure 2: Energy isosurfaces of the shell like solutions with K = l,m = 0fori? = l — 8. 
The images are coloured based on the value of 7:2 and the grey sphere represents the 
boundary of space in the Poincare ball model. 
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Figure 3: A plot of the energy for charge S = 1 — 5 shell like solutions against k 



kappa 


Figure 4: The numerical result of the energy compared to the rational map approximation 
for S = 2, for various value of k. If you consider the percentage of the approximation 
that the numerical result takes, it remains roughly constant within our numerical error. 











(a) cyrstal-chunk solution (b) shell-like solution 


Figure 5:5 = 8 static solution, (a) energy density plot of the crystal chunk solution with 
K = 1, m = 0, (b) energy density plot of the shell-like solution with k = 1 , m = 0 

4.2 Crystal chunk Solutions 

For the crystal chunk solutions we will consider a couple of cases, the 5 = 8 and 32 
solutions. In Euclidean space we find that the 5 = 8 solution needs a relatively high 
mass for the crystal chunk solution to become the global minima. This massive solution 
can be considered to be two 5 = 4 Skyrmions, joined along an axis perpendicular to 
a face of the shape. They have relative rotation of | around the axis joining the two 
solitons. 

In figure we observe that both the crystal chunk and shell-like solutions are attainable 
in hyperbolic space with k = 1. However, it appears that the crystal solution is the global 
minima for all non-zero curvatures considered. Note that the energies of the two solutions 
get very close and could be within numerical error of each other. The 5 = 8 crystal chunk 
solution is the lowest charge crystal solution and hence the energy difference might not 
be discernible with our accuracy. It is possible that the non-shell like solution does in fact 
become the minimal energy solution, for higher values of the curvature. We will consider 
a higher charge solution where the energy difference will be more discernible. Also, if 
the crystal chunk solutions act as with increasing the mass term in Euclidean space, we 
may find that the crystal chunk solution would become the minimal energy solution for 
a lower curvature. 

The 5 = 32 crystal chunk solution, displayed in figure |^b-c) , has a far lower energy 
than that of the shell like solution in figure |^a) for even low values of k. This is also the 
case for a small mass term in the Euclidean model. Hence we have demonstrated that 
not only are the profile functions related for Skyrmions with massive pions in Euclidean 
space and with massless pions in hyperbolic space, but the energetically favourable form 
of solution is also similar. 


5 Dynamics 


The solutions to the time-dependant equations of motion that follow from (2.1) were again 
found using a time dependent 4th-order Runga-Kutta method. The grid was modelled 
using the Poincare ball model of radius 1 (fixing k = 1) on a cubic grid with (201)^ grid 
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(a) S = 32 


(b) B = 32 


(c) B = 32 


Figure 6: Energy density plots of the multi-soliton solution for B = 32 for various 
isosurface values, coloured based on 112 value for (a) shell like solution with energy 40.43, 
(b-c) crystal chunk solution with energy 38.22. 



Figure 7: Scattering along a geodesic through the origin, with zero initial velocity, with 
solitons in the attractive channel (relative rotation of tt around a line perpendicular to 
the diagonal). 


point, hence the lattice spacing Ax = 0.005. Spatial derivatives have been approximated 
using a 4th-order finite difference method. The product ansatz was used for well separated 
single charge solitons. 

The simplest situation to consider is scattering along a geodesic that passes through 
the centre of the space, as seen in figure This gives a straight geodesic, with a clear 
parallel to Euclidean space and hence an obvious attractive channel (rotate relative by 
TT around an axis perpendicular to the connecting straight line). The scattering process 
then proceeds as expected with the solitons scattering at 7 r/ 2 . 

Due to hyperbolic translations (elements of the isometry group of hyperbolic space) one 
would expect in general, single Skyrmions to follow geodesics until they scatter. After 
scattering, the emerging Skyrmions will follow alternate geodesics, oriented to the incident 
paths by a rotation of tt around an orthogonal axis. The maximal channel will be a 
rotation of one of the solitons relative to the other by tt around an orthogonal axis to 
the tangent of the connecting geodesic. On scattering, the Skyrmions should merge to 
form the standard B = 2 solution, oriented to lie in the incident plane, however it may 
appear deformed due to the curvature of the space. The results presented here confirm 
these expectations and can be observed in figure 
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Figure 8 : Scattering along a curved geodesic, with zero initial velocity, in the attractive 
channel (relative rotation of vr around a line perpendicular to the geodesic). 

6 Conclusions 

We have found both static and dynamic solutions for hyperbolic Skyrmions of various 
curvature. The static solutions have been related to massive solutions in Euclidean space, 
by making use of the relation shown in [7j . It has been demonstrated that the link between 
curvature in hyperbolic space and mass in Euclidean space extends to full solutions of 
various topological charge, allowing predictions to be made for the type of solution that 
will occur in the two models. 

We have supplied evidence that suggests the rational map approximation is a good ap¬ 
proximation for increasing curvature. It seems to retain its accuracy regardless of the 
curvature considered. This would suggest that we can model Skyrmion solutions in the 
infinite curvature limit, by using their respective rational maps. This is analogous to the 
hyperbolic monopole case, where solutions for infinite curvature become rational maps 
[TT] . It would be interesting to see if there were some interesting limit in which it produces 
exact solutions, that in some way corresponds to hyperbolic monopoles. 

The dynamics of various soliton initial conditions have also been studied. The attractive 
channel was shown to be a relative rotation by vr around an axis orthogonal to the 
connecting geodesic. 

It would be interesting to consider the form of a soliton crystal in hyperbolic space, due 
to the interesting symmetries and tilings that can be formed from various polyhedron. 
It would be sensible to start with the 2-dimensional analogue, due to the difficulty of 
the task. Some similar work has been done with 2 -dimensional vortices in the hyperbolic 
plane, concentrating on the tiling with Schlafi symbol {8,8} [12]. 
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